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Abstract 

We study the impact of the four-quark dipenguin operator to the D° — D 
mixing. It is shown to contribute to the short distance piece at the same order 
of magnitude as the box diagram. 



1 



Typeset using REVTgX 



The phenomenon of meson mixing has been studied for a long time. Observed in the 
K° — K° and B° — B° systems, it provides an extremely sensitive test of the Standard Model 
(SM) and essential for the CP violation in the neutral meson system. 

It is well known that because of the GIM cancellation mechanism and large mass of the 
top quark the short distance box diagram dominates in B° — B° system and constitutes a 
significant fraction of the K°—K° mixing amplitude. The case of D°—D° system is somewhat 
special: the 6-quark contribution to the fermion loop of the box diagram providing AC = 2 
transition is diminished by a tiny V u b CKM matrix element. Thus, only the light quark 
mass difference guarantees that mixing does take place. The effect vanishes in the limit of 
the SU(3) invariance. All of that results in the estimated value for Amp being of the order 
of 10~~ 17 GeV if only short distance contributions are taken into account IJ, @. 

Calculating a box diagram and constructing the effective Hamiltonian one realizes that 
the smallness of the short distance piece is guaranteed by a factor of (m 2 — m^/M^m 2 [0: 

_G F a (m 2 - m 2 ) 2 

neff v/2 87rsin 2 ^ U il^m 2 
(«7 M (1 + 75)c«7^(l + 75)c + 2u(l- 7 5 )cu(l - 7 5 )c) (1) 

with & = V* c V iu . Here, the 6-quark contribution is dropped. This leads to the following 
expression for the Am^f [0], @ 

D x /2 47rsin 2 ^^ M 3 M^m 2 JD DK " D! ^0.2^ L/^J v > 

with ~ 132 MeV, fp ~ 165 MeV, and Bp = B' D = 1 in the usual vacuum saturation 
approximation to 

(D'lO^D ) = -£^°B D , (D°\0 2 \D°) = AV^Q>™k B < 
X 1 11 1 3 2m D ' 3 m 2 2m D D 

Ox = uj^l + 7 5 )cu7 M (l + 75) c, 2 = u(l - 75)cw(l - 7 5 )c (3) 

In contrast to the K-meson mixing, the appearance of the second operator can be traced to 
the fact that mass of the external c-quark provides a large momentum scale. As one can 
see from the Eq. (3.7) of |T|, the inclusion of the 6-quark further decreases the box diagram 
contribution. 



In this note we would like to address an additional contribution to the short- distance 
D-meson mixing amplitude which is topologically distinct from the box diagram - the so- 
called double penguin or "dipenguin" operator. This operator was initially introduced in |3| 
for K° — K° mixing amplitude and has been subsequently studied in Q in application to 
K° - K° as well as to B° - B° systems. It was shown to be marginally important in the 
former and completely negligible in the latter case. Here we will introduce the dipenguin 
operator for the AC = 2 transitions. It will be shown that this operator contributes to the 
D-meson mass difference at the same order of magnitude as the usual box diagram. 

The effective operator relevant to dipenguin AC = 2 transition can be obtained from 
the usual AC = 1 penguin vertex (we neglect a tiny dipole contribution): 

= J^ 7m (i + 75 )ycGr<9 2 - fr&)K (4) 

Here F\ is a modified Inami-Lim function ||. Using unitarity of the CKM matrix it reads 
Fi = J2i = £, S {F( — Ff) + C,b{Fi — Ff). It is common to discard 6-quark contribution to 
F\ as being suppressed by small V u b factors. Note, however, that by keeping it we enhance 
the Fx by approximately 20 — 30%. Also, for the intermediate 6-quark this vertex (as well as 
the following Hamiltonian) is local. From (Q) we obtain the following effective Hamiltonian 

H dp = -z^-F?({*Yv(l + i 5 )X a c)d^(u lu (l + 75 )A a c) - 

(u 7m (1 + 7 5)A a c)D(M7 M (l + 7 5 )A a c)) (5) 

In what follows, we denote two operators entering (|5|) as 0% and O2. In order to study the size 
of the dipenguin effects in charm mixing we derive an estimate of the Amo and compare 
with the usual box diagram contribution. To do that, in addition to the usual vacuum 
saturation approximation, we use pQCD in order to independently estimate the Am* by 
simply calculating the transition amplitude from Feynman diagram that determines the 
effect (Fig. la). We believe that pQCD can provide a reliable order-of-magnitude estimate of 
the dipenguin contribution since the momentum transferred through the gluon line in Fig. 
la is relatively large, Q 2 ~ — m? c . 
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Using the equations of motion and neglecting the up quark mass we obtain for the first 
operator in (|5|) 

Ox = U7„(l + 75 )A a c d»d v u lv {\ + 7s)A a c ~ m\ u{\ - 75 )A a c u{\ - 75 )A a c (6) 

Before we compare the estimate of the Am* with the relevant box diagram contribution 
we would like to note that the dipenguin diagram does not have a power dependence upon 
the internal quark masses. The leading behavior of the Inami-Lim function is logarithmic in 
m s ^), and the estimate of the operator brings about power dependence upon the external 
quark masses, i.e. m c . This feature distiguishes this contribution from that of the usual box 
diagram. Recalling the fact that the dominant contribution to K and B mixing amplitudes 
is proportional to the square of the top quark mass, it is not surprising that this effect is 
negligible in the K and B sectors. It is the fact of the "reduced" heavy quark dependence of 
the amplitude of D°—D° mixing which makes the dipenguin operator contribution effectively 
enhanced. 

Employing vacuum saturation method to estimate matrix elements we obtain 
(D^OAD ) = -&^-B, (D°\d 2 \D ) = 

v 1 11 1 9 2m D 9 2m D V ; 

with B and B' being the bag parameters. In addition to the vacuum saturation we assumed 
that each derivative acting on the quark field involves an average momentum of the quark. 
This yields 

Am* = 2(D°\H dp m = ^^F?(mlmim d )f D m D (ml - 4m c 2 ) (8) 

(2,71 71 

This formula deserves some additional discussion. In the case of D° — D° mixing the chief 

effect comes from the light quark sector. This is true for both the box and dipenguin 

diagrams, and makes the calculation a little more involved - one cannot simply discard the 

external quark momenta (masses). That is why one must use a modified expression for the 

Inami-Lim function || . The leading contribution to F{ comes from the integral 

r 1 r m 2 O 2 i rl m 2 O 2 i 

F{(mlQ 2 ) = -4 d x x{l - x)\n[-^ - ^x{l - x) = -4 - In + nft) (9) 

JO ^Vl w 1V1 W 1V1 W m i 
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where U(^) was defined in || for Q 2 > 0. In our case, of course, Q 2 < ~ — m 2 . The 

i 

necessity of the second term becomes obvious if one looks at the chiral limit, wherein the 
first logarithm in (^) blows up. The final result for F% involves a delicate cancellation among 
these contributions, yielding a final result \Fx\ = 0.01 — 0.02 (whereas Fx ~ 0(1) if the 
momentum flow is discarded). Comparing @ and @ we find 

I Am D I ^ a * F ?( m l m2 s, m d) m D (1Q) 
] Am b S xl ~ 8vr \Ud\(m 2 s - m 2 ) 2 1 ' 

where we have put m c m mo, and a s ~ 0.4. The relative size of the box and dipenguin 
contribution shows that the latter is of the same order of magnitude as the box diagram (Our 
estimate gives ~ 20 — 50% depending on the choice of quark masses). This is not surprising 
if one recalls that higher order QCD corrections tend to "smooth out" a power-like GIM 
suppression, just as in the case of 5-meson decays. This is not so relevant in the bottom (or 
strange) meson sector since a large mass of the top quark actually converts GIM-suppression 
to "GIM-enhancement" , thus making higher order corrections relatively unimportant. 

Another interesting observation is the fact that the dipenguin diagram actually con- 
tributes to the Am D with a sign opposite to the box diagram (compare (0) and (|8|))! This 
implies that the short- distance piece is even smaller than was claimed in previous estimates 
based solely on the box diagram contribution. Note that this operator gives rise to a whole 
family of the diagrams (e.g. Fig. lb) that by no means can be calculated in perturbative 
QCD but might be potentially important for D° — D° mixing. 

It must be stressed that D° — D° mixing is not dominated by the short distance box 
diagram contribution [[| but rather by long distance pieces 0, 0. This effect has been 
estimated in using dispersive techniques for a class of two-body pseudoscalar intermediate 
states and was shown to boost the value of Amp to Amp ~ 10~ 16 GeV. Indeed, it is not 
excluded that additional contributions could conspire in a way that they cancel the two 
body piece 0. 

In conclusion, we have estimated the contribution of the dipenguin diagram to the short- 
distance amplitude for D° — D° mixing. It is shown to contribute at the same order of 
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magnitude as the box diagram. 
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APPENDIX 

Here we will provide an independent pQCD estimate of the dipenguin matrix element 
(D°\O dp \D°). To do that we employ a Brodsky-Lepage exclusive QCD description || for 
calculating the Feynman diagram in Fig. la. In this formalism, the amplitude of interest can 
be expressed as a convolution of the mesons' distribution function with the hard scattering 
amplitude T(x,y). We assume that the momentum fraction carried by the c-quark in D 
meson is x, and the momentum fraction carried by the c quark in the D meson is 1 — y. As 
is well known, the heavy quark carries nearly all the momentum of the heavy-light bound 
state which makes the distribution functions strongly peaked at x ~ 1 and y ~ 0. This 
fact simplifies the choice of the form for the distribution amplitudes. In order to obtain the 
estimate of the effect we use the simplest form for the distribution amplitudes 

<Kz) = J2=5{l-x-e), </>{y) = J°=5{e-y) (11) 

These amplitudes are normalized such that 

[ <l>(x 1 )dx 1 = (0\A,\D°) = if D p Dlx (12) 

with xi — 1 — x. The desired mass difference is then given by 

Am% = 2(D°\H dp \D°) = -— f dxdy<f>* D (y)T(x,y)My) (13) 

771 £) J 

where a "scattering amplitude" T(x,y) can be read off the diagram shown in Fig. la and 
m£> in the denominator comes from the standard normalization of the meson states. A 
computation of Eq. (|13|) involves the effective vertex Vg M : 

V 8 , = ^Y8^( FliQ2) " Q "®\ (1 + 7s) " ™eW 2 )i<V^(l - 7 5 )) (14) 
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where Fi(Q 2 ) = J2k^kF ki (x k ,Q 2 ); k = b,s,d; x k = m 2 k /M^; F ki (x kl Q 2 ) are the correspond- 
ing modified Inami-Lim functions, and Q 2 is a momentum transferred through the gluon. 
Numerically F x = -0.015 for m b ~ 5 GeV, m s ~ 0.2 GeV, m d ~ 0.01 GeV. 
The calculation of (^) is relatively straightforward and yields 

lAmff^gm^^m-F,)"} (15) 

Here Cp = 4/3, a s ~ 0.4, and we put m c pa nip = 1.87 GeV. The calculation amounts to 
Am* ~ — 0.2- 10~ 17 (f D / f n ) 2 GeV . Note that the pQCD estimate and the vacuum saturation 
estimate give reasonably close results. 
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(a) (b) 
FIG. 1. Dipenguin diagram (a) and a possible long-distance contribution (b) 
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